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A novel theoretical approach to the dynamics analysis of excitation of plasmon modes in nanopar- 
ticles is presented. This approach is based on the biorthogonal plasmon mode expansion and it leads 
to the predictions of time-dynamics of excitation of specific plasmon modes as well as their steady 
state amplitude and their decay. Temporal characteristics of plasmon modes in nanoparticles are 
expressed in terms of their shapes, permittivity dispersion relations and excitation conditions. In the 
case of the Drude model, analytical expressions for time-dynamics of plasmon modes are obtained. 
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Plasmon resonances in nanoparticles have been the fo- 
cus of considerable experimental and theoretical research 
lately. This research is motivated by numerous scientific 
and technological applications of these resonances in such 
areas as near-field microscopy, nano-lithography, surface 
enhanced Raman scattering, nanophotonics, biosensors, 
optical data storage, etc. Steady state properties of plas- 
mon resonances in nanoparticles have been mostly stud- 
ied, while the dynamics of specific plasmon modes is the 
least understood area of plasmonics. This state of affairs 
has prompted the current burst of activity in the exper- 
imental research of plasmon dynamics. This research is 
exemplified by publications 0, 0, S II H, IS 0> 

where 

the plasmon dynamics and plasmon decay (dephasing) 
rates have been extensively studied by using advanced 
femtosecond techniques. The theoretical temporal anal- 
ysis of plasmon modes is also very important to fully 
comprehend the time-dynamics of mode excitation and 
decay as well as to estimate their steady state amplitude. 
The temporal analysis of plasmon resonance modes can 
be very instrumental in the area of light controllability of 
plasmon resonances in semiconductor nanoparticlesUQ, 
where proper time synchronizations of excitations of spe- 
cific plasmon modes may be needed. 

In this Letter, we present a self-consistent theoretical 
technique for the temporal analysis of specific plasmon 
modes in nanoparticles and compare (where it is possi- 
ble) the theoretical results with experimental data. This 
technique is based on the biorthogonal plasmon mode ex- 
pansion and it leads to the predictions of time-dynamics 
of excitation of specific plasmon modes and their decay as 
well as their steady state amplitude in terms of nanopar- 
ticle shapes, permittivity dispersion relations and excita- 
tion conditions. For instance, an explicit formula is de- 
rived for the steady state amplitude of plasmon modes in 
terms of real and imaginary parts of dielectric permittiv- 
ity, amplitude of incident field and its spatial orientation 
with respect to dipole moments of the plasmon modes. 



In the case of the Drude model, analytical expressions for 
plasmon dynamics are obtained which suggest that the 
reciprocal of the Drude damping factor can be identified 
with the decay (dephasing) time of plasmon modes. 

To start the discussion, consider a nanoparticle with 
boundary S and dielectric permittivity e(u>) in free-space 
with dielectric constant e . It has been demonstrated 
d, 0, US EH that the resonance values e k of dielectric 
permittivity and corresponding resonant plasmon modes 
can be found by solving the eigenvalue problems for the 
boundary integral equation 



a k (Q) = %L I ak (M) ^|^dS M) 
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where a k (M) has the physical meaning of surface elec- 
tric charges on S that produce electric field of k-th 
plasmon mode, T k (M) has the physical meaning of dipole 
densities on S that produce displacement field T) k of k-th 
plasmon mode, while all other notations have their usual 
meaning. The resonance values of e k and X k are related 
by 



£fc ~ £q 
£fe + £o ' 



(3) 



After X k are computed and resonance values of permit- 
tivity e k are found by using ([3]), resonance frequencies are 
determined from the dispersion relation: 



e fe = e'(uj k ) = Re [e(u) k )} . 



(4) 



It turns out that eigenfunctions a k {M) and t^(M) are 
biorthogonal: 



hi • 
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Thus, the set of eigenfunctions a k (M) can be used for 
the biorthogonal expansion of actual boundary charges 
a(M,t) induced on particle boundary during the excita- 
tion process: 



a(M,t) =J2 a k(t)<Tk(M), 



(6) 



fc=i 



where, according to (O, the expansion coefficient ak(t) is 
given by the formula: 



a k (t)=j> a(M,t)Tk(M)dS M - 



(7) 



It is clear that the time evolution of the expansion coef- 
ficient dfe(t) reveals the time-dynamics of fc-th plasmon 
mode corresponding to the eigenfunction a k (M). Since 
the medium of resonant nanoparticle is dispersive and ex- 
hibits nonlocal in time constitutive relation D(t) vs E(f), 
the frequency domain technique will be employed for the 
calculations of ctk(t). This means that the equation ©- 
([7]) will be Fourier transformed 



er(M,w) = ^a k {u>)a k {M), 



fc=l 



a k (iv) 



a(M,co)r k (M)dS M , 
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and the equation for a k (u>) will be first derived and 
solved. Subsequently, a k (t) will be found through inverse 
Fourier transform. 

To derive the equation for d k (uj), the Fourier trans- 
formed boundary condition for the normal components 
of electric field on S is invoked: 

e(w)E+(Q,uj)-e E-(Q,uj) = [e -e(uj)] E^{Q,u). (10) 

Here: E+(Q : uj) and E~(Q,ui) are the limiting values of 
the normal components of Fourier transformed electric 
field at Q € S from inside and outside S, respectively, 
while En^ (Q, to) is the Fourier transformed normal com- 
ponent of the incident field on S. 

By using formula a(Q,uj) = e [E-(Q,u)) - E+(Q,uj)}, 
the boundary condition (fTU|) can be rearranged as follows: 

cr(Q,w) 



e(w) - e 
Now it can be recalled [l 
a(Q,w) 



E+(Q,u) = EW(u> 



(11) 



hj that 



E+(Q,w) = - 
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By substituting (Tj"2")) into pT|) and then by using ex- 
pansion ([5]) and formulas JTJ and ©, the following ex- 
pression for a k (oS) is derived: 



The last integral can be simplified, because the incident 
field E( )(Q,t) is practically uniform within nanoparti- 
cles. Thus, this field can be represented as E^°'(Q,t) — 
E /(<). By using this form of E(°)(Q, t) in UHJ) and tak- 
ing into account formula and the fact that due to the 
"scaling" freedom provided by normalization condition 
(j5J the following expression § s r\QT k (Q)dSQ = (ek-e )v k 
for the dipole moment pfc of k-th plasmon mode (see [9j) 
can be used, the last formula can be represented in the 
form: 



where 



a k {u) = (E • p fe )3fcH/(w), 
e(w) - e 



e k - e(w) ' 



(14) 
(15) 



and f(uj) is the Fourier transform of f(t). 

Formulas (fT4|) -(fl~5 |) are remarkably simple and they 
clearly suggest that g k (co) can be construed as normalized 
(by Eo - pfc) transfer function of fc-th plasmon mode. The 
formula (|15p reveals the resonance nature of excitation of 
fc-th plasmon mode at the frequency u) k . Indeed, accord- 
ing to at this frequency e k — e(u k ) — —ilm [e(ui k )] and 
the magnitude of g k (u>) will be narrow peaked if e"{u> k ) 
is sufficiently small, that is when the specific plasmon 
resonance is strongly pronounced. 

By using formula pT]) . the following expression is de- 
rived for a k (t): 

a k (t) = (Eo • p*) / g k (t-t')f(t')dt', (16) 
Jo 

with g k (t) being the inverse Fourier transform of g k {Lo). 
Thus, the algorithm of analysis of time-dynamics of k- 
th plasmon mode can be stated as follows. First, the 
eigenvalue problem ([I]) and ([3]) is solved and resonance 
values e k of dielectric permittivity and surface electric 
charges a k (M) of the corresponding plasmon mode are 
found. Next, by using a k (M), the dipole moment p,t is 
computed. Then, g k (t) is determined through the inverse 
Fourier transform of g k (w) given by (|14[) . Finally, formula 
(fTBf is employed to evaluate the time evolution of a k (t) 
which reveals the time-dynamics of fc-th plasmon mode. 
It is worthwhile to mention that the outlined computa- 
tions can be performed by using actual, experimentally 
measured e(u>). 

Expressions (fT4|l - (fT6|) can be useful for analytical cal- 
culations as well. To demonstrate this, consider the 
steady state case when f(t) — sinuj k t, where uj k is the 
resonance frequency defined by the formula Then, 
f(u) = i^f\ [5(uj - uj k ) - 8(uo + uj k )] and, by using ([14])- 
(|15[) as well as the inverse Fourier transform of a k (u) , we 
arrive at: 



d k (cj) 



2e Q X k [e(uj) - epj 
2e A fe + [e(w)-e ] (A fc -1) 



E$\Q, u;)T k (Q)dS Q . a[ ss \t) - -(E • p k ) 
(13) 



e'(oj k ) - e 
e"(u k ) 



cosuj k t + sinuj k t 



(17) 
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where superscript (ss) indicates the steady state of afc(t). 

In the case of strong plasmon resonances when \e'(ujk) — 
eo| >> |e"(a>fc)|, the last formula is simplified as follows: 



Ass) 



(t) 



-(E • Pfe 



COSLUkt- 



(18) 



As typical for most resonances, the steady state is 
shifted by almost ninety degree in time with respect to 
the incident field Eo(i). It is also natural that the magni- 
tude of the steady state is controlled by the ratio of real 
and imaginary parts of dielectric permittivity at the res- 
onance frequency. It is also revealing that the resonance 
magnitude depends on the spatial orientation of the in- 
cident field Eo(f) with respect to the dipole moment pfe 
of plasmon resonance mode. 
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FIG. 1: C(wo) computed for: (1) gold nanocylinders placed 
on a glass substrate that resonate at wavelength 774 nm (see 
i) and (2) gold nanorings placed on a glass substrate that 
resonate at wavelength 1102 nm (see [3). The Au dispersion 
relation from [l5| have been used in computations. 



For off-resonance excitation f(t) = sintoot, similar cal- 
culations lead to the expression: 



a i SS) W = ( E o ' Pk)C(uj )cos(uj t + tp), 



where 



C(w ) 



' [e'M-e ] 2 + [e"MP 
[e k - e'M] 2 + [e"M] 2 



(19) 



(20) 



Figure 1 illustrates that C(ujq) is narrow peaked at the 
resonance frequency u>k for the two cases of plasmon res- 
onances experimentally studied in Q and [3] . 

Next, we point out that simple analytical expressions 
for g k (t) in (fT6)) can be obtained for the Drude model of 



e(w) = e 



^7, 



U){ljJ + £7) 



(21) 
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FIG. 2: Dynamics of plasmon resonances for Au nanorings on 
a glass substrate computed for the resonance wavelength 1102 
nm (see 0) by using (a) Durde model (7 = 1.0753 x 10 14 s" 1 
[l5l ]) and (b) Au dispersion relation from [TBI ]. 



(22) 



Indeed, by using (l21j) in (115]) we obtain: 



g k (u) = - 



[iu> — otx)(iui — a.2) ' 



where: ol\ 



-i/3 and a 2 = %-i0, (3 = V 4 "^ 2 



Now 



the inverse Fourier transform of (|22p yields for t > 



9k(t) 



+ 7 2 ) e - ¥ 



sin/3t. 



(23) 



Formulas (fTrJ)) and (f2"3"|) can be used for analytical cal- 
culations of a,k(t) for "rectangular" laser pulses E(°)/(t). 
Indeed, in the case when: 



/(*)) 



for t < 0, 



(24) 



1 sinWfci for t > 0, 
from (T6J), and (23) we derive: 

o fc (t) = -(E -p fe )w fc e-** ficos/3t-| Sl n/5^+4 ss) (t), 

(25) 
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where 



■"' s) (t) = (E • Pfe) ( —cosuj k t - sinuj k t ) 



(26) 
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In the case of finite rectangular laser pulse 



/(*) 



from OB, J23J) and 




for t < 0, 

for T < t < 0, 

for t > T, 



we derive for t > T: 



(27) 



a k (t) = -(E • Pfc)e ^ 



-i(t-T) 



uj k e 



— cosuj^t — sinuj^t 

1 



1 2 

—cos/3t — —sin/3t 

1 P 

(28) 



It is instructive to compare the computational results 
obtained by using the Drude model with those obtained 
by using the experimentally measured dispersion rela- 
tion e(u>). This comparison is presented by Figures 2a 
and 2b for Au rings subject to finite rectangular laser 
pulses. These figures suggest that the Drude model leads 
to quantitatively similar results as the use of actual dis- 
persion relation. Formula ([28]) also implies that I/7 can 
be identified as the decay (dephasing) time for the light 
intensity of plasmon modes. In accordance with the avail- 
able data for 7 (see 0, EH), this suggests that formula 
(f2"8"| predicts the decay (dephasing) time for plasmon 
modes in gold (and silver) nanoparticles in the range of 
5-12fs, which is consistent with the experimental results 
reported in 0, 0, |, Q, 0, |, 0] . 



mon mode (thin line). The computed plasmon time- 
dynamics was used to compute the third order autocor- 
relation function (ACF) which was compared with the 
experimentally measured ACF from Q. The results of 
comparison presented in Figure 4. Figure 5 presents the 
comparison between our calculations of the second-order 
ACF for noncentrosymmetric L-shape gold nanoparticles 
from Q and the experimentally measured in [|| second- 
order ACF. Figures 4 and 5 suggest the agreement with 
experimental data within ten percent. 
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FIG. 4: Envelop of the calculated third order ACF (filled cir- 
cles) superimpose on the measured third order ACF (solid 
line)from 0] for Au cylinders on a glass substrate at wave- 
length 774 nm (maximum is normalized to 32). 
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FIG. 3: Dynamics of plasmon resonances for Au cylinders on 
a glass substrate computed for the resonance wavelength 774 
nm (see 0]). The Au dispersion relation from [l5| have been 
used in computations. 

Finally, we present the comparison between computa- 
tional results based on the technique outlined above and 
experimental results presented in @, 0| . By using equa- 
tions p5]) - pB we have computed the time-dynamics of 
the 774 nm resonance wavelength plasmon mode for Au 
cylinders on glass substrate [2|. Figure 3 presents the 
time variation of the incident electric field of the laser 
pulse (bold line) used in experiments reported in and 
the corresponding computed time-dynamics of the aver- 
age (over nanoparticle volume) electric field of the plas- 
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FIG. 5: The calculated second order ACF superimpose on 
the envelop of measured second order ACF from [(J for non- 
centrosymmetric L-shape Au nanoparticles on a tantalum- 
dioxide substrate at wavelength 838 nm (maximum is nor- 
malized to 8). 
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